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METHOD FOR LOAD BALANCING A LOOP OF 
PARALLEL PROCESSING ELEMENTS 



BACKGROUND OF THE INVENTION 

[0002] The present invention relates generally to parallel processing and more particularly to 

balancing the work loads of the processing elements within a parallel processing system. 

[0003] Conventional central processing units ("CPU's"), such as those found in most 

personal computers, execute a single program (or instruction stream) and operate on a single 

stream of data. For example, the CPU fetches its program and data from a random access 

memory ("RAM"), manipulates the data in accordance with the program instructions, and 

writes the results back sequentially. There is a single stream of instructions and a single 

stream of data (note: a single operation may operate on more than one data item, as in X = Y 

+ Z, however, only a single stream of results is produced). Although the CPU may determine 

the sequence of instructions executed in the program itself, only one operation can be 

completed at a time. Because conventional CPUs execute a single program (or instruction 

stream) and operate on a single stream of data, conventional CPUs may be referred to as a 

single-instruction, single data CPU or an SISD CPU. 
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[0004] The speed of conventional CPUs has dramatically increased in recent years. 
Additionally, the use of cache memories enables conventional CPUs faster access to the 
desired instruction and data streams. However because conventional CPUs can complete only 
one operation at a time, conventional CPUs are not suitable for extremely demanding 
applications having large data sets (such as moving image processing, high quality speech 
recognition, and analytical modeling applications, among others). 

[0005] Improved performance over conventional SISD CPUs may be achieved by building 
systems which exhibit parallel processing capability. Typically, parallel processing systems 
use multiple processing units or processing elements to simultaneously perform one or more 
tasks on one or more data streams. For example in one class of parallel processing system, 
the results of an operation from a first CPU are passed to a second CPU for additional 
processing, and from the second CPU to another CPU, and so on. Such a system, commonly 
known as a "pipeline", is referred to as a multiple-instruction, single-data or MISD system 
because each CPU receives a different instruction stream while operating on a single data 
stream. Improved performance may also be obtained by using a system which contains many 
autonomous processors, each running its own program (even if the program running on the 
processors is the same code) and producing multiple data streams. Systems in this class are 
referred to as a multiple-instruction, multiple-data or MIMD system. 
[0006] Additionally, improved performance may be obtained using a system which has 
multiple identical processing units each performing the same operations at once on different 
data streams. The processing units may be under the control of a single sequencer running a 
single program. Systems in this class are referred to as a single-instruction, multiple data or 
SIMD system. When the number of processing units in this type of system is very large (e.g., 
hundreds or thousands), the system may be referred to as a massively parallel SIMD system. 
[0007] Nearly all computer systems now exhibit some aspect of one or more of these types 
of parallelism. For example, MMX extensions are SIMD; multiple processors (graphics 
processors, etc) are MIMD; pipelining (especially in graphics accelerators) is MISD. 
Furthermore, techniques such as out of order execution and multiple execution units have 
been used to introduce parallelism within conventional CPUs as well. 
[0008] Parallel processing is also used in active memory applications. An active memory 
refers to a memory device having a processing resource distributed throughout the memory 
structure. The processing resource is most often partitioned into many similar processing 
elements (PEs) and is typically a highly parallel computer system. By distributing the 
processing resource throughout the memory system, an active memory is able to exploit the 
very high data bandwidths available inside a memory system. Another advantage of active 
memory is that data can be processed "on-chip" without the need to transmit the data across i 
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system bus to the CPU or other system resource. Thus, the work load of the CPU may be 
reduced to operating system tasks, such as scheduling processes and allocating system 
resources. 

[0009] A typical active memory includes a number of interconnected PEs which are capable 
of simultaneously executing instructions sent from a central sequencer or control unit. The 
PEs may be connected in a variety of different arrangements depending on the design 
requirements for the active memory. For example, PEs may be arranged in hypercubes, 
butterfly networks, one-dimensional strings/loops, and two-dimensional meshes, among 
others. 

[0010] In typical active memories, load imbalances often occur such that some PEs are idle 
(i.e., without assigned tasks) while other PEs have multiple tasks assigned. To maximize the 
effectiveness of the active memory, it is desirable to balance the work load across all of the 
PEs. For example in an active memory having a multitude of identical PEs, it is desirable that 
each PE be assigned the same number of instructions by the central sequencer, thus 
maximizing the resources of the active memory. Additionally in an active memory having 
non-identical PEs, it may be desirable to assign more tasks to the PEs with greater processing 
capabilities. By balancing the load, the amount of time that one or more PEs is idle while 
waiting for one or more other' PEs to complete their assigned tasks is minimized. 
[0011] Thus, there exists a need for a method for balancing the load of a parallel processing • 
system such that the resources of the parallel processing system are maximized. More 
specifically, there exists a need for a method for balancing the load of an active memory such 
that the resources of the active memory are maximized. 

SUMMARY OF THE INVENTION 

[0012] One aspect of the present invention relates to a method for balancing the load of a 
parallel processing system having a plurality of parallel processing elements arranged in a 
loop, wherein each processing element has a local number of tasks (v r ) associated therewith, 
wherein r represents the number for a selected processing element PE r , and wherein each of 
the processing elements are operable to communicate with a clockwise adjacent processing 
element and with an anti-clockwise adjacent processing element. The method comprises 
determining within each processing element (PE r ) a total number of tasks (V) present within 
the loop, calculating a local mean number of tasks (M r ) within each of the plurality of 
processing elements (PE r ), and calculating a local deviation (D r ) within each of the plurality 
of processing elements (PE r ). The method also comprises determining a sum deviation within 
each of the processing elements (PE r ) for one -half the loop in an anti-clockwise direction (A) 
wherein the one-half of the loop being relative to each of the selected processing elements 
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(PE r ), determining a sum deviation within each of the processing elements (PE r ) for one-half 
the loop in a clockwise direction (Q wherein the one-half of the loop being relative to each of 
the selected processing elements (PE r ), determining a clockwise transfer parameter (T c ) and 
an anti-clockwise transfer parameter (T a ) within each of the processing elements (PE r ), and 
redistributing tasks among the plurality of processing elements in response to the clockwise 
transfer parameter (T c ) and the anti-clockwise parameter (T a ) within each of the plurality of 
processing elements (PE r ). 

[0013] Another aspect of the present invention relates to a method for assigning tasks among 
a plurality of processing elements within a parallel processing system, the processing 
elements being connected in a loop and having a local number of tasks associated therewith. 
The method comprises determining the total number of tasks on the loop, computing a local 
mean value for each of the processing elements, and computing a local deviation for each of 
the processing elements, the local deviation being representative of the difference between the 
local number of tasks for a processing element and the local mean value for the processing 
element. The method also includes summing the deviation of the processing elements located 
within one-half of the loop in an anti-clockwise direction, summing the deviation of the 
processing elements located within one-half of the loop in a clockwise direction, computing a 
number of task to transfer in a clockwise direction, computing a number of task to transfer in 
an anti-clockwise direction, and redistributing tasks relative to the number of task to transfer 
in a clockwise direction and the number of task to transfer in an anti-clockwise direction. 
[0014] The present invention enables tasks to be distributed along a group of serially 
connected PEs so that each PE typically has X number of tasks or (X+l) number of tasks to 
perform in the next phase. The present invention may be performed using the hardware and 
software (i.e., the local processing capability) of each PE within the array. Those advantages 
and benefits, and others, will become apparent from description of the invention below. 

BRIEF DESCRIPTION OF THE DRAWINGS 

j 

[0015] To enable the present invention to be easily understood and readily practiced, the 
present invention will now be described for purposes of illustration and not limitation, in 
connection with the following figures wherein: 

[0016] FIG. 1 is a block diagram illustrating an active memory according to an embodiment 
of the present invention. 

[0017] FIG. 2 is a block diagram of a processing element for the active memory illustrated in 
FIG. 1 according to an embodiment of the present invention. 

[0018] FIG. 3 illustrates an array of the processing elements illustrated in FIG. 2 arranged in 
a loop according to an embodiment of the present invention. 
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[0019] FIG. 4 illustrates an operational process for balancing the load within a loop of 
processing elements according to various embodiments of the present invention. 
[0020] FIG. 5 illustrates the determination of the sum deviation in the anti -clockwise half of 
loop (A) and the sum deviation in the clockwise half of loop (Q is determined for a local PE 
according to an embodiment of the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

[0021] As discussed above, parallel processing systems may be placed within one or more 
classifications (e.g., MISD, MIMD, SIMD, etc.). For simplicity, the present invention is 
discussed in the context of a SIMD parallel processing system. More specifically, the present 
invention is discussed in the context of a SIMD active memory. It should be noted that such 
discussion is for clarity only and is not intended to the limit the scope of the present invention 
in any way. The present invention may be used for other types and classifications of parallel 
processing systems. 

[0022] FIG. 1 is a block diagram illustrating an active memory 10 according to an 
embodiment of the present invention. It should be noted that the active memory 10 is only 
one example of a device on which the methods of the present invention may be practiced and 
those of ordinary skill in the art will recognize that the block diagram of FIG. 1 is an 
overview of an active, memory device 10 with a number of components known in the art being 
omitted for purposes of clarity. 

[0023] Active memory 10 is intended to be one component in a computer system. 
Processing within active memory 10 is initiated when the active memory 10 receives 
commands from a host processor (not shown), such as the computer system's CPU. A 
complete processing operation (i.e., data movement and processing) in the active memory 10 
may consist of a sequence of many commands from the host to the active memory device 10. 
[0024] Active memory 10 is comprised of a host memory interface ("HMI") 12, a bus 
interface 14, a clock generator 16, a task dispatch unit ("TDU") 18, a DRAM control unit 
("DCU") 20, a DRAM module 22, a programmable SRAM 24, an array control sequencer 26, 
and a processing element array 28, among others. 

[0025] The HMI 12 provides an input/output channel between the host (such as a CPU, not 
shown) and the DRAM module 22. In the current embodiment, the HMI 12 receives 
command (cmd), address (addr), and data signals (among others) from and sends data and 
ready (rdy) signals (among others) to the host. The HMI 12 approximates the operation of a 
standard non-active memory so that the host, without modifications, is compatible with the 
active memory 10. 
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[0026] The HMI 12 may be similar in its operation to the interface of a synchronous DRAM 
as is know in the art. Accordingly, the host must first activate a page of data to access data 
within a DRAM module 22. In the current embodiment, each page may contain 1024 bytes of 
data and there may be 16,384 pages in all. Once a page has been activated, it can be written 
and read through the HMI 12. The data in the DRAM module 22 may be updated when the 
page is deactivated. The HMI 12 also sends control signals (among others) to the DCU 20 
and to the processing element array 28 via the task dispatch unit 18. 

[0027] The HMI 12 may operate at a frequency different than that of the frequency of the 
master clock. For example, a 2x internal clock signal from clock generator 16 may be used. 
Unlike a traditional DRAM, the access time for the HMI 12 uses a variable number of cycles 
to complete an internal operation, such as an activate or deactivate. Thus the ready signal 
(rdy) is provided to allow the host to detect when a specific command has been completed. 
[0028] The bus interface 14 provides and input/output channel between the host and the 
TDU 18. For example, the bus interface 14 receives column select (cs), write command (w), 
read command (r), address (addr), and data signals (among others) from and places interrupt 
(intr), flag, and data signals (among others) onto the system bus (not shown). The bus 
interface 14 also receives signals from and sends signals to TDU 18. 

[0029] The clock generator 16 is operable to receive an external master clock signal (xl) and 
operable to provide the master clock signal (xl) and one or more internal clock signals (x2, 
x4, x8) to the components of the active memory. It should be apparent to one skilled in the 
art that other internal clock signals may be produced by the clock generator 16. 
[0030] The TDU 18 communicates with the bus interface 14, the HMI 12, the programmable 
SRAM 24, the array control sequencer 26, and the DCU 20. In the current embodiment, the 
TDU 1 8 functions as an interface to allow the host to issue a sequence of commands to the 
array control sequencer 26 and the DCU 20. Task commands from the host may be buffered 
in the TDU's FIFO buffers to allow a burst command to be issued. Commands may contain 
information on how the tasks in the array control sequencer 26 and the DCU 20 should be 
synchronized with one another, among others. 

[0031] The DCU 20 arbitrates between the TDU 18 and the HMI 12 and sends commands to 
the DRAM modules 22 and the processing element array 28. The DCU 20 also schedules 
refreshes within the DRAM modules 22. In one embodiment, the DRAM modules 22 of the 
active memory 10 may be comprised of sixteen 64k xl28 eDRAM (or embedded DRAM) 
cores. Each eDRAM core may be connected to an array of sixteen. PEs, thus providing 256 
(16 x 16) PEs in all. 

[0032] The programmable SRAM 24 functions as a program memory by storing commands 
issued by the TDU 18. For example, the TDU 18 may transmit a "write program memory 
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address" command which sets up a start address for a write operation and a "write program 
memory data" command which writes a memory location and increments the program 
memory write address, among others. The programmable SRAM 24, in the current 
embodiment, has both an address register and a data output register. 
[0033] The array control sequencer 26 may be comprised of a simple 16 bit minimal 
instruction set computer (16-MISC). The array control sequencer 26 communicates with the 
TDU 1 8, the programmable SRAM 24, and the DCU 20, and is operable to generate register 
file addresses for the processing element array 28 and operable to sequence the array 
commands, among others. 

[0034] The processing element array 28 is comprised of a multitude of processing elements 
("PEs") 30 (see FIG. 2) connected in a variety of different arrangements depending on the 
design requirements for the processing system. For example, processing units may be 
arranged in hypercubes, butterfly networks, one-dimensional strings/loops, and two- 
dimensional meshes, among others. In the current embodiment, the processing elements 30 
are arranged in a loop (for example, see FIG. 3). The processing element array 28 
communicates with the DRAM module 22 and executes commands received from the 
programmable SRAM 24, the array control sequencer 26, the DCU 20, and the HMI 12. Each 
PE in the processing element array 28 includes dedicated H-registers for communication with 
the HMI 12. Control of the H-registers is shared by the HMI 12 and the DCU 20. 
[0035] Referring now to FIG. 2, a block diagram of a PE 30 according to one embodiment of 
the present invention is illustrated. PE 30 includes an arithmetic logic unit ("ALU") 32, Q- 
registers 34, M-registers 36, a shift control and condition register 38 (also called "condition 
logic" 38), a result register pipeline 40, and register file 42. The PE 30 may also contain other 
components such as multiplexers 48 and logic gates (not shown), among others. 
[0036] In the current embodiment, the Q-registers 34 are operable to merge data into a 
floating point format and the M-Registers 36 are operable to de-merge data from a floating 
point format into a single magnitude plus exponent format. The ALU 32 is a multiplier-adder 
operable (among others) to receive information from the Q-registers 34 and M-registers 36, 
execute tasks assigned by the TDU 18 (see FIG. 1), and transmit results to the shift control 
and condition logic 38 and to the result register pipeline 40. The result register pipeline 40 is 
operable to communicate with the register file 42, which holds data for transfer into or out of 
the DRAM modules 22 via a DRAM interface 44. Data is transferred between the PE and the 
DRAM module 22 via a pair a registers, one register being responsive to the DCU 20 and the 
other register being responsive to the PE 30. The DRAM interface receives command 
information from the DCU 20. The DRAM interface 44 also permits the PE 30 to 
communicate with the host through the host memory access port 46. 
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[0037] In the current embodiment, the H-registers 42 are comprised of synchronous SRAM 
and each processing element within the processing element array 28 contains eight H-registers 
42 so that two pages can be stored from different DRAM locations, thus allowing the 
interleaving of short i/o bursts to be more efficient. Result register pipeline 40 is also 
connected to one or more neighborhood connection registers ("X-register") (not shown). The 
X-register links one PE 30 to its neighboring PE's 30 in the processing element array 28. 
[0038] The reader desiring more information about the hardware shown in FIGs. 1 and 2 is 
directed to UK Patent application (serial no. not yet assigned) entitled "Control of Processing 
Elements in Parallel Processors" filed 17 September 2002, (Micron no. 02-1604) which is 
hereby incorporated by reference. Details about the PEs may also be found in UK Patent 
Application No. 021562.2 entitled "Host Memory Interface for a Parallel Processor" filed 17 
September 2002, (Micron no. 02-0703) which is hereby incorporated by reference. 
[0039] FIG. 3 is a simplified diagram showing the interconnections of an array of the 
processing elements 30 (as illustrated in FIG. 2) arranged in a loop 50 according to an 
embodiment of the present invention. In the current embodiment, loop 50 is comprised of 
eight (8) PEs 30 (i.e., PEo, PEi, . . . PE 7 ) which are interconnected via their associated X- 
register links. It should be noted that the number of PEs 30 included in loop 50 may be 
altered while remaining within the scope of the present invention. As illustrated in FIG. 3, 
each PE is operable to communicate with its clockwise, and anti-clockwise neighbor. For 
example, PEi is operable to communicate with its clockwise neighbor, PE 2 , and with its anti- 
clockwise neighbor, PE 0 . In the current embodiment, every PE 30 on the loop 50 receives 
instructions from a single TDU 18 as discussed in conjunction with FIG. 1. Furthermore, 
each PE has a local number of tasks (v r ) associated therewith. For example, PEo has three (3) 
tasks associated therewith (i.e., v 0 = 3), PE t has six (6) tasks associated therewith (i.e., v ; = 6), 
PE 2 has two (2) tasks associated therewith (i.e., v 2 = 2), etc. 

[0040] FIG. 4 illustrates an operational process 60 for balancing the work loads between the 
PEs 30 on loop 50 according to an embodiment of the present invention. Operational process 
60 begins by determining the total number of tasks (V) present on the loop in operation 61. 
As discussed above in conjunction with FIG. 3, each PE r (where r = 0 to 7, e.g., PE 0j PEi, . . . 
PE 7 ) in the loop has a local number of tasks (v r ) associated therewith. In the current 
embodiment, each PE r passes its own value v r onto its clockwise neighbor and simultaneously 
receives a value v r+/ from its anti -clockwise neighbor. Each PE r keeps a running partial sum 
(i.e., adds each value v r+/ received to its own value v r ). This process continues until each 
value v r has moved clockwise around the loop and visited each PE r , in this case seven 
transfers are needed. At the end of the rotation process, the sum represents the total number 
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of tasks (V) on the loop. As illustrated in FIG. 3, loop 50 has forty-three (43) total tasks 
associated therewith. 

i=N-\ 

[0041] The sum (V) can be expressed by the equation V = £ v y , where //represents the 

/=o 

number of PEs 30 in the loop 50 (here N = 8), and v, represents the local number of tasks 
associated with an i th processing element in the loop. For example, for i = 3, the number of 
tasks associated with PE 3 (i.e., v 3 ) is added to the sum V. It should be noted that after a 
rotation is completed, each PE r will have calculated the same value for (V). It should also be 
noted in the current discussion, "local" refers to the values or functions associated with a 
single PE within the loop, whereas "global" refers to the values or function associated with 
the entire loop of PEs. 

[0042] After the total number of tasks (V) present on the loop is determined in operation 61, 
the local mean number (M r ) of tasks for each PE r is computed in operation 62. In the current 
embodiment, operation 62 employs a rounding function to ensure that no tasks are lost or 

/"=/V-I 

"gained" during the rounding process (i.e., to ensure that V = ). The rounding 

(=0 

function is discussed in more detail in U.S. Patent Application Serial No. entitled 

"Method for Rounding Values for a Plurality of Parallel Processing Elements" filed 

(DBOO 1064-000, Micron no. 02-1269) and incorporated in its entirety by reference herein. 
[0043] For example assume that 13 tasks (i.e., V= 13) are to be shared by the eight PEs (i.e., 
PE 0 through PE 7 ). Without the rounding function, the local mean for each PE would be PE r = 
1 .625 before rounding (i.e., 13- 8 = 1 .625). If the fraction thirteen-eighths is set to round 
down for each PE (i.e., 13 - 8 = 1), then the sum of the means for all of the individual PEs 
(i.e., PE 0 through PE 7 ) is equal to eight (8) and five (13 -8 = 5) tasks are lost. In contrast, if 
the fraction thirteen-eighths is set to round up for each PE (i.e., 13-8 = 2), then the sum of 
the means for all of the individual PEs (i.e., PE 0 through PE 7 ) is equal to sixteen (16) and 
three (16-13 = 3) extra tasks are gained. 

[0044] The rounding function M r = Trunc((V + E r )/N) prevents tasks from being lost or 
gained (where M r represents the local mean for PE r , W represents the total number of PEs 30 
in the loop 50, and E r represents a number in the range of 0 to (N-\)). In the current 
embodiment, each PE is assigned a different E r value for controlling the rounding. The 
simplest form for the function E is the case in which E r = P n where P r represents the PEs 
position round the loop. For example, for PE 0 , E 0 = 0; for PE,, E, = 1; for PE 2 , E 2 = 2; etc. 
By assigning each PE in the loop a different E r value, the rounding function can be controlled 
such that some of the local means are rounded up and some of the local means are rounded 
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down, thus insuring that V = / M Mj • It should be noted that in the current embodiment, the 

1=0 

local mean for each PE 30 in the loop is computed in parallel with the local means of the 
other PEs in the loop. 

[0045] Table 1 illustrates the local mean calculation for the loop 50 as illustrated in FIG. 3 in 
which the total number of tasks on the loop is equal to forty -three (43). Referring to Table 1, 
it is apparent that the rounding function controls the rounding such that M 0 through M 4 are all 
rounded to five (5), whereas M 5 through M 7 are all rounded to six (6) . The sum of the values 
of M 0 through M 7 is equal to forty-three (43), which equals the total number of tasks (V) on 
the loop. Thus, tasks are neither lost nor gained due to rounding. 
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(V+E r )/N 
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Table #1 - Local Mean Calculation for the Loop 50 (V= 43, N = 8). 

I 

[0046] After the local means are computed in operation 62, the local deviation D r is 
calculated for each PE in operation 63. In the current embodiment, the local deviation is 
simply the difference between the local number of tasks and the local mean (i.e., D r = v r - 
M r ). The local deviations for PE 0 through PE 7 are illustrated in Table #1 . 
[0047] After the local deviations are computed in operation 63, the sum deviation in the anti- 
clockwise half of loop (A) is determined for each PE in operation 64. The anti-clockwise sum 
04) is found in a similar manner as that used to determine the partial value sum (V) in 
operation 61. In operation 64, however, the local deviations (D r ) are only rotated halfway 
around the loop in clockwise direction and summed. The sum deviation in the anti-clockwise 

half of loop can be represented by the equation: A = ^Z),- . 

i=i 

[0048] After the sum deviation in the anti-clockwise half of loop (^4) is determined in 
operation 64, the sum deviation in the clockwise half of loop (C) are determined for each PE 



0O433599.DOC 



10 



in operation 65. The clockwise sum (Q is found in a similar manner as that used to 
determine the anti-clockwise sum (A) in operation 64. In operation 65, however, the local 
deviations (D r ) are rotated halfway around the loop in an anti-clockwise direction and 
summed. The sum deviation in the clockwise half of loop can be represented by the equation: 

i=(N /2)+l 

[0049] FIG. 5 illustrates how the sum deviation in the anti-clockwise half of loop (A) and the 
sum deviation in the clockwise half of loop (Q is determined for PE 2 . As seen in FIG. 5, the 
sum deviation in the clockwise half of loop (Q may be illustrated as a combination of PE i3 
PE,, and PE 5 into a "super PE". Referring to Table #1, the sum deviation of this super PE is 
thus D 3 + D 4 + D 5 = 2 + 3 + (-1) = 4 (i.e., C = 4). Likewise, the sum deviation in the anti- 
clockwise half of loop (A) may be illustrated as a combination of PE y , PE 0 , and PE 7 into 
another "super PE". Referring to Table #1, the sum deviation of this super PE is thus D, + D 0 
+ £> 7 = 1 +(-2) +1=0 (i.e., ,4 = 0). 

[0050] After the sum deviation in the clockwise half of loop (Q is determined in operation 
65, clockwise and anti-clockwise transfer parameters {T c and T a , respectively) are determined 
in operation 66. In the current embodiment, T c is determined using the equation T c = (25- C 
+ A)/4 and T a is determined using the equation T a = (2S-A + Q/4. 

[0051] Referring again to FIG. 5 from the perspective of PE 2 , the loop has four values C, A, 
0, and S, where C represents both the sum deviation in the clockwise half of loop and the 
deviation of the first "super PE", A represents both the sum deviation in the anti-clockwise 
half of loop and the deviation of the second "super PE", S represents the deviation of the 
selected PE (e.g., here PE 2 ) and O represents the deviation of the PE opposite to the selected 
PE (i.e., the PE for which the local deviation is being determined; here, PE/). The selected PE 
can deduce the deviation value of its opposite PE) because all deviations in the loop must sum 
to zero (i.e., ,4 + C+ S+ O = 0, thus O = -^-C-5'=-0-4 - (-3) « -1). It should be noted 
that A and C are calculated for each PE in parallel. Thus in the current embodiment, the 
clockwise and anti-clockwise transfer parameters for PE 2 are T c = [2(-3) - 4 + 0]/4 = -2.5 and 
r fl = [2(-3)-0 + 4)/4 = -0.5. 

[0052] It should be noted that the values obtained for T c and T a may need to be rounded in 
such a manner that R(T C ) + R{T a ) = D r . In the current embodiment, tasks are transmitted in 
only one direction at a time around the loop (i.e., either in the clockwise or anti-clockwise 
direction). A direction is selected for the 'first' transmission around the loop and the values 
for T c and T a are rounded up in this direction. It should be noted that by ensuring 'excess 
traffic' is sent in the 'first' direction, the chance of the process finishing one step earlier is 
increased. In the current embodiment, tasks are transmitted in the anti-clockwise first, such 
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that R(T a ) = Ceil (T a ), where the 'Ceil ' function returns the closest integer greater than or 
equal to the supplied input. To ensure that extra tasks are not created or lost by the rounding 
of R(T a l R(T C ) is set equal D r - R(T a ). 

[0053] Accordingly in the example above, tasks are transmitted anti-clockwise first such that 
T a = -0.5 is rounded up to 0. To ensure that extra tasks are not created or lost by the rounding 
of R(T a ), R(T C ) is set equal D r - R(T a ). Thus, T c is equal to -3 (i.e., -3- 0 = -3). It should be 
noted that other rounding mechanisms may be used while remaining within the scope of the 
present invention. For example, T c may be rounded up on odd numbered PE's and T a rounded 
up on even numbered PE's such that pairs of odd and even PE's exchange their 'excess 
traffic'. * 
[0054] In the case where the loop 50 is comprised of an odd number of PEs 30, an extra 
"phantom" PE may be used. The phantom PE is assigned a deviation of zero (0) and is 
located diametrically opposite from the perspective of the selected PE. For example, assume 
that PE 7 is removed from loop 50 such that only seven PEs remain (i.e., PE 0 to PE 6 ). To 
calculate the local deviation of PE 0 , the phantom PE would be placed between PE 3 and PE 4 ; 
for PE l3 between PE 4 and PE 5 ; for PE 2 , between PE 5 and PE 6 , etc. Thus, the phantom PE is 
located in the loop 50 such that the number of PEs between the selected PE and the phantom 
PE in the clockwise direction is equal to the number of PEs between the selected PE and the 
phantom PE in the anti-clockwise direction. 

[0055] In some instances, it may be desirable to clamp the transfer rates to reduce the 
number of iterations needed to balance the loop. In one embodiment, a non-linear clamping 
operation is utilized. For example, the equations for T c and T a may be re-written as T c = 
Trunc[(2S + A) - 4] and T a = Trunc[(2S -; A) - 4], respectively, where A = (A-C) represents 
the number of 'thru' tasks (i.e., the number of tasks passing through the current PE. If T c and 
T a are of opposite sign, then the number of 'thru' values may be reduced by clamping either 
T c or T a to zero. The remaining value (i.e., T a or T c respectively) may then be found using the 
identity S = T a + T c . This ensures that any rounding error introduced by the Trunc function is 
correctly compensated for such that S is finally equal to zero. 

[0056] In the current embodiment for example, the transfer parameter with the smallest 
absolute magnitude may be selected. The desired result can be achieved by applying the 
following non-linear modifications to A where Mag = abs(2S). If A > Mag then set A equal to 
Mag and if A < -Mag, then set A equal to -Mag. The revised value for A (i.e., A = Mag or A 
= -Mag) is then substituted into the equations T c = Trunc[(2S + A) - 4] and T a = Trunc[(2S - 
A) - 4] . It should be noted that other clamping operations may be used while remaining 
within the scope of the present invention. 
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[0057] After the clockwise and anti-clockwise transfer parameters are determined in 
operation 66, the tasks are redistributed among the PEs in response to the clockwise and anti- 
clockwise transfer parameters (i.e., T c and T a , respectively) in operation 67. In the current 
embodiment, a positive T c parameter represents the number of values that are to be 
transmitted clockwise out of the local PE. A negative T c parameter represents the number of 
values that are to be transmitted from the clockwise PE into the local PE. Similarly, a 
positive T a parameter represents the number of values that are to be transmitted anti- 
clockwise out of the local PE. A negative f a parameter represents the number of values that 
are to be transmitted from the anti-clockwise PE into the local PE. 
[0058] If the local deviation (D) is negative, one or more of the received values will be 
"absorbed" by the local PE to make up the local deficit. The others will be transmitted, either 
from the clockwise PE to the anti-clockwise PE, or from the anti-clockwise PE to the 
clockwise PE. On occasion, some PEs may start off with no values at all, these PEs may have 
to "mark time" until they receive a value. It should be noted that after each successful 
transmission or receipt, the local parameters T c and T a need to be updated. The redistribution 
stage only terminates when T c = T a = 0 for all PEs. 

[0059] As discussed above, the clockwise transfer parameter for PE 2 in the current 
embodiment is T c = -3. Because T c is negative, 3 tasks are to be transmitted from PE 3 into 
PE 2 . Likewise, the anti-clockwise transfer parameter for PE 2 in ihe current embodiment is T a 
= 0. Because T a is zero, no tasks are to be transmitted from PE 7 into PE 2 . It should be 
apparent that PE 2 had a deviation of -3 (i.e., D 2 = -3). Thus, three tasks were transferred into 
PE 2 in operation 67. 

[0060] It should be recognized that the above -described embodiments of the invention are 
intended to be illustrative only. Numerous alternative embodiments may be devised by those 
skilled in the art without departing from the scope of the following claims. 
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What is claimed is: 

1 . A method for balancing the load of a parallel processing system having a plurality of 
parallel processing elements arranged in a loop, wherein each processing element (PE r ) has a 
local number of tasks associated therewith, wherein r represents the number for a selected 
processing element, and wherein each of said processing elements is operable to communicate 
with a clockwise adjacent processing element and with an anti -clockwise adjacent processing 
element, the method comprising: 

determining within each processing element a total number of tasks present 
within said loop; 

calculating a local mean number of tasks within each of said plurality of 
processing elements; 

calculating a local deviation within each of said plurality of processing elements; 

determining a sum deviation within each of said processing elements for one-half 
said loop in an anti-clockwise direction, said one -half of said loop being relative to 
each of said selected processing elements; 

determining a sum deviation within each of said processing elements for one-half 
said loop in a clockwise direction, said one-half of said loop being relative to each of 
said selected processing elements; 

determining a clockwise transfer parameter and an anti-clockwise transfer 
parameter within each of said processing elements; and 

redistributing tasks among said plurality of processing elements in response to 

said clockwise transfer parameters and said anti-clockwise parameters within each of 

s 

said plurality of processing elemjents. 

I 

2. The method of claim 1 wherein said determining within each of said processing 
elements a total number of tasks present within said loop, comprises: 

transmitting said local number of tasks associated with each of said processing 
elements to each other of said plurality of processing elements within said loop; 

receiving within each of said processing elements said number of local tasks 
associated with said each other of said plurality of processing elements; and 

summing said number of local tasks associated with each of said processing 
elements with said number of local tasks associated with each other of said plurality 
of processing elements. 
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3. The method of claim 1 wherein said determining said total number of tasks present 

/=iV-I 

within said loop includes solving the equation V = ^]v, , where N represents the number of 

processing elements in said loop, and v,- represents said local number of tasks associated with 
an i th processing element in said loop. 

4. The method of claim 1 wherein said calculating a local mean number of tasks within 
each of said plurality of processing elements includes solving the equation 

M r = Trunc{{V + E r )/N), where M r is said local mean for PE r , N is the total number of 
processing elements in said loop, and E r is a number in the range of 0 to (Af-l) and wherein 
each processing element has a different E r value. 

5. The method of claim 3 wherein E r controls said Trunc function such that said total 
number of tasks for said loop is equal to the sum of the local mean number of tasks for each 

of said plurality of processing elements in said loop (i.e., V = ^M,. ). 

/=o 

6. The method of claim 3 wherein said local mean M r = Trunc((V + E r )/ N) for each 
local PE r within said loop is equal to one of and (AM-1). 

7. The method of claim 1 wherein said calculating a local deviation within each of said 
plurality of processing elements, comprises finding the difference between said local number 
of tasks and said local mean number for each of said plurality of processing elements. 

8. The method of claim 1 wherein said determining a sum deviation within each of said 
processing elements for one -half of said loop in an anti -clockwise direction comprises: 

transmitting said local deviation associated with each of said processing elements 
halfway around said loop in an anti -clockwise direction, said one-half of said loop 
being relative to each of said selected processing elements; 

receiving said local deviation associated with each other of said processing 
elements halfway around said loop in a clockwise direction, said one-half of said 
loop being relative to each of said selected processing elements; and 

summing said local deviations associated with each other of said processing 
elements halfway around said loop in a clockwise direction. 
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9. The method of claim 1 wherein said determining a sum deviation within each of said 
processing elements in one-half of said loop in a clockwise direction comprises: 

transmitting said local deviation associated with each of said processing elements 
halfway around said loop in an clockwise direction, said one-half of said loop being 
relative to each of said selected processing elements; 

receiving said local deviation associated with each other of said processing 
elements halfway around said loop in a anti-clockwise direction, said one-half of said 
loop being relative to each of said selected processing elements; and 

summing said local deviations associated with each other of said processing 
• elements half way around said loop in an anti-clockwise direction. 

10. The method of claim 1 wherein said determining a clockwise transfer parameter and 
an anti-clockwise transfer parameter within each of said processing elements comprises: 
setting the clockwise transfer parameter equal to (2S + A - Q + 4; and 
setting the anti-clockwise transfer parameter equal to (25 + C - A) + 4, where S 
represents the deviation of a selected processing element; C represents said sum 
deviation in a clockwise half of loop relative to said selected processing element, and 
A represents said sum deviation in an anti-clockwise half of loop relative to said 
selected processing element. 

11. The method of claim 1 wherein said determining a clockwise transfer parameter and 
an anti-clockwise transfer parameter within each of said processing elements comprises at 
least one of: 

setting the clockwise transfer parameter equal to Trunc[(2S + A) - 4] and setting 
the anti-clockwise transfer parameter equal to S - T c and 

setting the anti-clockwise transfer parameter equal to Trunc[(2S - A) + 4] and 
setting the clockwise transfer parameter equal to S - T a ; 
where A = Mag, if A > Mag, where A = -Mag, if A < -Mag, where Mag = abs(2S), and where 
S represents the local deviation of a selected processing element. 

12. A method for reassigning tasks among an odd numbered plurality of processing 
elements within a parallel processing system, said processing elements being connected in a 
loop and each having a local number of tasks associated therewith, the method comprising: 

determining the total number of tasks on said loop; 

computing a local mean value for a selected processing element; 
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computing a local deviation for said selected processing element, said local 
deviation representative of the difference between said local number of tasks for said 
selected processing element and said local mean value for said selected processing 
element; 

inserting a phantom processing element within said loop; 

summing said deviation of said processing elements located within one-half of 
the loop in an anti-clockwise direction relative to said selected processing element; 

summing said deviation of said processing elements located within one-half of 
the loop in a clockwise direction relative to said selected processing element; 

computing a number of tasks to transfer in a clockwise direction for said selected 
processing element; 

computing a number of tasks to transfer in an anti-clockwise direction for said 
selected processing element; and 

reassigning tasks relative to the said number of tasks to transfer in a clockwise 
direction and said number of task to transfer in an anti-clockwise direction. 

13. The method of claim 12 wherein said determining the total number of tasks on said 
loop, comprises: 

transmitting said local number of tasks associated with each of said processing 
elements to each other of said plurality of processing elements within said loop; 

receiving within each of said processing elements said number of local tasks 
associated with said each other of said plurality of processing elements; and 

summing said number of local tasks associated with each of said processing 
elements with said number of local tasks associated with each other of said plurality 
of processing elements. 

14. The method of claim 12 wherein computing a local mean value for a selected 
processing element includes solving the equation M r = Trunc((V + E r )l N), where M r is 
said local mean for a PE r , N\s the total number of processing elements in said loop, and E r is 
a number in the range of 0 to (N-\). 

15. The method of claim 14 wherein E r controls said Trunc function such that said total 
number of tasks for said loop is equal to the sum of the local mean number of tasks for each 
of said plurality of processing elements in said loop and wherein each processing element has 
a different E r value assigned. 
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1 6. The method of claim 1 2 wherein said inserting a phantom processing element within 
said loop further comprises: 

locating said phantom processing element in a position within said loop that is 
diametrically opposed to said processing element; and 

assigning a zero deviation value to said phantom processing element. 

17. The method of claim 12 wherein said computing a local mean value for a selected 
processing element, said computing a local deviation for said selected processing element, 
said inserting a phantom processing element within said loop, said summing said deviation of 
said processing elements located within one-half of the loop in an anti-clockwise direction, 
summing said deviation of said processing elements located within one -half of the loop in a 
clockwise direction, computing a number of tasks to transfer in a clockwise direction for said 
selected processing element, computing a number of tasks to transfer in an anti-clockwise 
direction for said selected processing element, and reassigning tasks relative to the said 
number of task to transfer in a clockwise direction and said number of tasks to transfer in an 
anti-clockwise direction are completed simultaneously for each of said plurality of processing 
elements within said loop. 

18. The method of claim 12 wherein said summing said deviation of said processing 
elements located within one-half of the loop in an anti-clockwise direction relative to said 
selected processing element comprises: 

transmitting said local deviation associated with each of said processing elements 
half way around said loop in an anti-clockwise direction, said one-half of said loop 
being relative to each of said selected processing elements; 

receiving said local deviation associated with each other of said processing 
elements half way around said loop in a clockwise direction, said one-half of said 
loop being relative to each of said selected processing elements; and 

summing said local deviations associated with each other of said processing 
elements half way around said loop in a clockwise direction. 

19. The method of claim 12 wherein summing said deviation of said processing elements 
located within one-half of the loop in a clockwise direction relative to said selected processing 
element comprises: 

transmitting said local deviation associated with each of said processing elements 
half way around said loop in an clockwise direction, said one-half of said loop being 
relative to each of said selected processing elements; 
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receiving said local deviation associated with each other of said processing 
elements halfway around said loop in a anti-clockwise direction, said one-half of said 
loop being relative to each of said selected processing elements; and 

summing said local deviations associated with each other of said processing 
elements halfway around said loop in an anti-clockwise direction. 

20. A memory device carrying a set of instructions which, when executed, perform a 
method comprising: 

determining within each processing element a total number of tasks present 
within said loop; 

calculating a local mean number of tasks within each of said plurality of 
processing elements; 

calculating a local deviation within each of said plurality of processing elements; 

determining a sum deviation within each of said processing elements for one-half 
said loop in an anti-clockwise direction, said one -half of said loop being relative to 
each of said selected processing elements; 

determining a sum deviation within each of said processing elements for one -half 
said loop in a clockwise direction, said one-half of said loop being relative to each of 
said selected processing elements; 

determining a clockwise transfer parameter and an anti-clockwise transfer 
parameter within each of said processing elements; and 

redistributing tasks among said plurality of processing elements in response to 
said clockwise transfer parameters and said anti-clockwise parameters within each of 
said plurality of processing elements. 
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ABSTRACT OF THE DISCLOSURE 
A method for balancing the load of a parallel processing system having a plurality of 
parallel processing elements arranged in a loop, wherein each processing element has a local 
number of tasks associated therewith comprising determining within each processing element 
a total number of tasks present within the loop, calculating a local mean number of tasks 
within each processing element, and calculating a local deviation within each processing 
element. The method also comprises determining the sum deviations within each processing 
element for one-half the loop in an anti-clockwise direction and in a clockwise direction, 
determining clockwise and anti-clockwise transfer parameters within each processing 
element, and redistributing tasks among the processing elements in response to the clockwise 
and anti-clockwise transfer parameters. 
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